Towards unambiguous calculation of the topological entropy for mixed states 
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Calculation of topological order parameters, such as the topological entropy and topological mu- 
tual information, are used to determine whether states possess topological order. Their calculation 
is expected to give reliable results when the ground states of gapped Hamiltonians are considered, 
since non-topological correlations are suppressed by a finite correlation length. However, studies of 
thermal states and the effects of incoherent errors require calculations involving mixed states. Here 
we show that such mixed states can effectively lead to a diverging correlation length, and hence 
may give misleading results when these order parameters are calculated. To solve this problem, 
we propose a novel method to calculate the quantity, allowing topologically ordered states to be 
identified with greater confidence. 

PACS numbers: 03.65.Ud, 03.67.Mn, 71.10.Pm, 73.43.Nq 
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Topologically ordered states have generated a great 
deal of interest, both in condensed matter physics [lj- 
5] and quantum information theory. This is due in the 
most part to their ability to support anyonic quasipar- 
ticles, useful for fault-tolerant quantum computation [Q- 
03. Topologically ordered states cannot, in general, be 
detected by any local order parameter. Instead order pa- 
rameters such as the topological entropy and topological 
mutual information must be used [ 1 lM 1 31 ] . These aim to 
detect the unique non-local correlations found in topolog- 
ically ordered states. However, present formulations are 
known to also have a contribution from non-topological 
correlations [III 03 • When pure ground states of gapped 
Hamiltonians are considered, these non-topological corre- 
lations only extend for a short range determined by their 
correlation length. Their contributions are therefore ex- 
pected to be suppressed when the regions used to evalu- 
ate the entropies are made larger than this length scale. 
However, here we show that statistical mixtures of the 
degenerate ground states of gapped Hamiltonians with 
local interactions can effectively have a diverging correla- 
tion length, preventing the suppression of non-topological 
contributions. This could cause non-topological states to 
be identified as topological, or topological states to be 
identified as non-topological. This is an important mat- 
ter to address, since calculations using mixed states are 
required for the studies of thermal states and incoherent 
error models that are of high physical and practical sig- 
nificance [HI, Il5l - fl7| . To solve the problem, a modified 
definition of the topological entropy is proposed. This 
correctly identifies and discards the non-topological con- 
tributions, and so can be used with confidence the de- 
termine whether topological order is present for mixed 
states. 

To begin, it will be advantageous to introduce the con- 
cepts which will be used throughout this work. The set 
of anyonic quasiparticles supported by a given topologi- 
cally ordered state can be described in terms of an anyon 
model. Within each of these models the quantum dimen- 
sion d a can be defined for each species of anyon, a jf| 
The total quantum dimension of the model is then given 



by D 2 = ^ a d\. The case of D = 1 corresponds to the 
trivial anyon model, which contains only the vacuum sec- 
tor, and is associated with non-topological states. The 
case of D > 1, however, is a signature of topologically 
ordered states. 

The realization of topological order requires large 
many-body systems. We use p to denote the state of 
such a system and pr = tr# e p to denote the state of 
the subsystem contained within an arbitrary region R, 
where R c is its complement. Using the Von-Neumann en- 
tropies of the states of these regions (for example Sr — 
— tr(p^lnpfl)) the mutual information between two re- 
gions R and R' may be defined as Zr,h' = Sr+Sr>— Srri 
[H ] . Here RR' is used as shorthand for the composite re- 
gion RU R' . We will use the convention that, if only a 
single region is specified in the subscript, the mutual in- 
formation is between that region and it's complement, 
e.g. I R = Ir.r c . 



I. THE TOPOLOGICAL ORDER PARAMETERS 

The two main methods to calculate a topological order 
parameter for mixed states are the topological entropy 
of Levin and Wen [121 ] (and a variant by Castelnovo and 
Chamon [n| ) and the topological mutual information of 
Iblisdir et al. [HI (based on the pure state topologi- 
cal entanglement entropy of Kitaev and Preskill [111]). 
Both of these measures result in the same value of 2 log D 
for topologically ordered states associated with an anyon 
model of total quantum dimension D. 

The method of Levin and Wen notes that states of 
topologically ordered systems have correlations for re- 
gions that form closed loops that do not exist for those 
along open loops 03 • As such an entropic difference be- 
tween the two cases is expected. This can be measured 
by taking an annulus and dividing it into two regions, A 
and B, then further dividing B into B\ and B2, as shown 
in Fig. [Ha). Using these, 
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(Sb — Sbx) — {Sab — Sab ± )- 
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Here the ellipsis represents contributions from non- 
topological correlations, which will be suppressed as the 
regions are made larger than the correlation length. In 
this work, we note that this calculation of the topolog- 
ical entropy may alternatively be expressed in terms of 
mutual informations as, 
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Note that this is not a modification of the quantity, just 
an alternative way of expressing it. By expressing it in 
this way we see that it is the information that one sec- 
tion of a closed string (A) shares with the rest (B), but 
does not share with its neighbouring sections (Bi). This 
is because A is correlated to B by both topological and 
local correlations, but is correlated to B\ by only local 
correlations. The difference between these mutual infor- 
mations thus isolates the topological correlations. There 
is no reason to suppose that this argument applies only 
to pure states. The entropy as defined by Levin and Wen 
is therefore equally applicable to mixed states. As such, 
this entropy has been applied to mixed states in previous 
studies [la - [l7j . It should be noted that a modified defini- 
tion for mixed states has also been proposed in [15| . but 
it gives equivalent results and so will not be considered 
explicitly here. 

The topological mutual information of Iblisdir et al. is 
a generalization of the topological entanglement entropy 
of Kitaev and Preskill [11] to mixed states. It uses the 
fact that the mutual information between the region R 
and its complement, i? c , will satisfy an area law of the 
following form, 



Ir = aL 



li 



Here a is a positive constant and the correction 7 is due 
to topological correlations. The ellipsis represents contri- 
butions from non-topological correlations, which disap- 
pear as the regions are made larger than the correlation 
length. 

To cancel out the boundary terms and isolate the topo- 
logical term, the following linear combination of mutual 
informations is taken using the regions shown in Fig. 
IHa), 

7/ + • • • = ~Ia — Ib — Ic — I ABC 

+ Iab+Iac + Ibc- (3) 

Note that the quantity called the topological mutual in- 
formation in [? ] is actually —7/. However, 7/ is used 
here such that the values of both measures considered 
are positive for topologically ordered states. With this 
definition, topologically ordered states associated with 
an anyon model of total quantum dimension D will have 
Ilw = li — 2 log 2, and non-topologically ordered states 
will have jlw = 7/ = 0. 

It is important to note that jlw and 7/ are measures 
of topological correlations, not of entanglement. For pure 
states, the distinction is unnecessary, since the presence 
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FIG. 1. The regions used for the calculation of a: (a) the 
topological entropy in Eq. @; (b) the topological mutual 
information in Eq. (3). For both of these entropy calculations, 
the size of the regions should be taken to be large for the most 
accurate results. 



of topological (or any other) correlations implies that en- 
tanglement is present. However, for mixed states, topo- 
logical correlations may arise from classical probability 
distributions, as noted previously in [l^]. This property 
can cause some ambiguity, since it is possible to confuse 
'true' topological order which is quantum in nature with 
classical topological order. However, it is not this prob- 
lem that we address here. 



II. LONG-RANGE NON-TOPOLOGICAL 
CORRELATIONS 

The corrections present in the calculations above are 
due to non-topological correlations that may be present 
in the system. These can be characterized by by two- 
point correlation functions of the form Gij = (OiOj) — 
(Oi){Oj), where Oi is an operator defined in the neigh- 
bourhood of a point i of the system. It is known 
that, for (pure) ground states of gapped Hamiltonians, 
Gi.j — 0(e~ d ^), where d is the distance between points i 
and j. As such, these correlations are suppressed beyond 
the finite correlation length £. Making all length scales 
for the regions used in the above calculations much larger 
than £ therefore successfully suppresses the corrections. 
However, difficulties can arise when mixed states are con- 
sidered. Though Gi.j decays with distance for each in- 
dividual degenerate ground state of a gapped Hamilto- 
nian, it does not necessarily do so for a mixture of them. 
The mixed state can therefore, effectively, have an infi- 
nite correlation length. Such long-range non-topological 
correlations will then cause the corrections to remain fi- 
nite for regions of arbitrary size, and the results of the 
calculations in Eq. and Eq. ([3]) to give an ambiguous 
characterization of topological order. 

As we will show below, this effect is not restricted to 
unrealistic states. Instead, it can occur for ground states 
of local gapped Hamiltonians. Furthermore, since the 
equally weighted mixture of degenerate ground states de- 
scribes the thermal state at zero temperature, it is an 
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FIG. 2. The spin lattice on which the toric code, and hence 
the states p' 1 ' and p^ 2 \ is defined. Dots denote spins, and 
operators are defined around each plaquette and vertex. The 
lattice has periodic boundary conditions. 



important and realistic state to consider for any Hamil- 
tonian. As such, it is vitally important to determine a 
definition of the topological entropy that is immune to 
such effects. 

Two simple examples of states with long-range non- 
topological correlations can be defined using the toric 
code model Q with a modified Hamiltonian. This model 
uses the spin lattice of Fig. [5J defining operators A v — 
Yii£v a f 011 t ne spins around each vertex v, and operators 
Bp = Ylitp °t around the spins of each plaquette p. The 
usual Hamiltonian of the model is, 

v p 

The A v and B p operators define the occupation of 
anyons. States within their —I (+1) eigenspace denote 
that an anyon is present (not present) on the correspond- 
ing vertex or plaquette. The Hamiltonian assigns energy 
to anyons, and hence the ground states of the Hamil- 
tonian are those with no anyons anywhere. Any eigen- 
state of the Hamiltonian is topologically ordered, with 
Ilw = 7i — 2 log 2. The same is true for any mixture 
of eigenstates if they all correspond to the same anyon 
configuration. As such, the ground states of the Hamil- 
tonian have Jlw = li = 2 log 2, as does any mixture of 
them. We will use p^ to denote their equally weighted 
mixture. 

For mixtures of Hq eigenstates with different anyon 
configurations, the topological order parameters can take 
any value from to 2 In 2, inclusive. The value Jlw — 
7/ = 2 log 2 is only obtained if the following condition 
is met [lg. Consider an annulus, such as that formed 
by the region AB in Fig. Q] (a), jlw = 7/ = 2 log 2 
if, by local operations on the spins within the annulus, 
a state can be prepared such that the parity of anyons 
on both the plaquettes and vertices enclosed by the the 
annulus (i.e. those with either full or partial support on 
spins within the region enclosed by the annulus) is in a 
definite state. 

Let us now consider the following modified version of 



the toric code Hamiltonian, 

Hi = — J A v A v i — J BpB p /. 

(v.v') <p,p') 

Here (v, v') denotes a pair of vertices joined by an edge, 
and (p,p') denotes a pair of plaquettes which share an 
edge. Hence, rather than single vertex and single plaque- 
tte terms of Hq, the Hamiltonian is composed of near- 
est neighbour two- vertex and two-plaquette interactions. 
Nevertheless, it is important to note that it is still both 
gapped and local, and has the same eigenstates as Hq. 
The ground states of Hi correspond to the states for 
which the either no vertices hold an anyon, or all vertices 
hold an anyon. The same is true for the plaquettes. Let 
us consider the equally weighted mixture of all ground 
states, p^ 1 '. 

In this case, the natural choice for the operator Oi for 
the correlation function Gij is A v (alternatively, B p ). 
It can then be easily seen that G v y(p( l >) = 1 for any 
pair of vertices v and v', without any decay due to the 
distance between them. The state therefore effectively 
has an infinite correlation length, despite being composed 
of ground states of a gapped and local Hamiltonian. 

It is clear that this state is topologically ordered, with 
Ilw =7/ = 2 log 2, since measurement of even a single 
vertex and plaquette on an annulus is sufficient to deter- 
mine the occupancies of all vertices and plaquettes. The 
parity of both vertex and plaquette occupations enclosed 
by the annulus is then both definite and known. 

To perform the calculation of the topological mutual 
information, note that state p^ can be decomposed as 
an equally weighted mixture of four mixed states. Each 
of these corresponds to a different anyon configuration 
(vertices and plaquettes all empty, vertices filled and pla- 
quettes empty, vertices empty and plaquettes filled, and 
vertices and plaquettes all filled), and hence they are or- 
thogonal. Also, the entropies for each of these mixed 
states are equal, since they are equivalent by local uni- 
taries. The reduced density matrices p^ and p^ are 
similarly composed of four orthogonal and locally uni- 
tarily equivalent mixed states, since the different anyon 
configurations can also be distinguished and manipulated 
within each region. The behaviour of the Von Neumann 
entropy when applied to block diagonal matrices [l3[ then 
means that the mutual information can be decomposed 
into two contributions. The first comes from the correla- 
tions within each of the orthogonal and locally equivalent 
states. Since these states are all equivalent to p(°\ this 
contribution is Jr(/jW). The second contribution comes 
from the fact that there are four states, distinguishable 
in both regions, that are equally mixed. This adds an 
amount 2 log 2 to the total mutual information. Hence, 

/K(p (1) ) = /fl,fl c (p (0) ) + 21og2. (4) 

When these mutual informations are combined accord- 
ing to Eq. ([3]), the contribution of In(p^) from each 
will lead to a total contribution of 2 log 2. This is due to 
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the topological correlations within the state p^°\ How- 
ever, the total contribution of the 2 log 2 terms will be 
—2 log 2. In total, therefore, the calculation for the topo- 
logical mutual information yields 7/ + . . . = 0. 

This result may lead anyone calculating this quantity 
to believe that p^ does not possess topological corre- 
lations, and therefore that the nearest neighbour inter- 
actions of H\ are not sufficient to support topologically 
ordered states. However, such a conclusion would be in- 
correct. The result of 7 + . . . = arises in this case only 
because non-topological correlations are not suppressed, 
no matter how large the regions are made, and are there- 
fore able to cause ambiguity. 

In general, any state on a 2D surface for which the spins 
(or other physical medium) around one point have cor- 
relations with other those around points spread through- 
out the surface, and the correlations do not decay with 
distance, will cause ambiguities in the calculation of the 
topological mutual information. If pure states occur with 
the same properties, similar ambiguities will arise for the 
topological entanglement entropy that it is based upon 
[111 ] . This is because such correlations cause the mutual 
information across any bipartition to have contributions 
that do not depend on the length of the boundary, and 
hence cannot be distinguished from genuine topological 
contributions. 

For the second example we consider another modified 
toric code Hamiltonian, 

H2 = — J A V A V > — J BpBp>. 

(v,v') (p,p') 

Here (v, v') denotes a pair of vertices joined by a verti- 
cal edge, and (p,p') denotes a pair of plaquettes which 
share a horizontal edge. As such, H2 is like Hi in that 
it consists of nearest neighbour 2- vertex and 2-plaquette 
interactions. But, unlike H±, these interactions are con- 
fined to a specific direction. The ground state of H 2 con- 
sist of any state for which each vertical line of vertices or 
plaquettes is either empty or full of anyons. We consider 
the equally weighted state p^ of all ground states. 

As before, using Oi = A v yields G VjV >(pW) — 1 for any 
pair of vertices v and v' on the same line. There is again 
no decay due to distance, and hence an infinite correla- 
tion length despite being composed of ground states of 
the a gapped and local Hamiltonian. 

This state is again topologically ordered with jlw — 
7/ = 2 log 2. This is because measurement of a vertex 
(plaquette) on each line of vertices (plaquettes) that in- 
tersect the region enclosed by the annulus is sufficient to 
determine the parity of the vertex (plaquette) occupa- 
tions within this region. 

Like pd\ the state p^ is a sum of orthogonal and lo- 
cally equivalent states of different anyon configurations, 
with each state locally equivalent to p(°K The same is 
true for the reduced states for any region. The entropy 
of a region R is therefore composed of an Sr(p^) con- 
tribution due to the correlations within each state and a 
contribution of log 2 for each vertical line of plaquettes 



or vertices that intersects the region (i.e. has at least one 
vertex or plaquette with full support on the region). If 
we use n R to denote the number of such lines, the mutual 
information between two regions is 

Ir,R' 

When these are combined according to Eq.([2]), the result 
will be jlw + ■ ■ • = 2 log 2 + nlog2. Here n denotes 
the number of vertical lines that intersect both A and 
B2 but do not intersect Bi . The value of the topological 
entropy is therefore overestimated by an amount that is 
not suppressed as the size of the regions are increased, 
but in fact increases. This is due to the non-topological 
correlations along lines of anyons. The same is true for 
the reformulation of the Levin and Wen method in [l5| . 

In general, any states with infinite range non- 
topological correlations along lines are able to cause am- 
biguities to the topological entropy of Levin and Wen. 
This is because such correlations are able to directly cor- 
relate the regions A and B2, and the calculation is built 
on the assumption that such correlations do not occur. 

With these examples, we see that both of these topo- 
logical order parameters are susceptible to long range 
non-topological correlations, since they can cause correc- 
tions that cannot be suppressed. This may then lead to 
non-topological phases being identified as topological by 
mistake, or vice versa, especially as more exotic systems 
are probed (20l - [23 |. Though both examples considered 
only fool one of the parameters, with p^> and p^ both 
giving the expected values of 2 In 2 when using jlw & n d 
7/ respectively, clearly states will exist which fool them 
both. As such, a way to rid the entropy of the non- 
topological corrections must be determined. 

III. MODIFYING THE TOPOLOGICAL 
ENTROPY 

To solve the problem, let us take a closer look at how 
the topological entropy is defined. In Eq. © it was 
shown that jlw is the information that a section of a 
closed string (^4) shares with the rest of the string (B), 
but not with the neighbouring sections (Bi). Since it 
can usually be expected that any contributions to Ia,b 
due to local correlations act over the boundary between 
A and its neighbour Bi, this definition serves to remove 
all local boundary correlations and thus, in most cases, 
measures only the topological correlations around the 
string. The problem, as mentioned above, is that suf- 
ficiently long range local correlations contribute in a way 
that does not depend on the boundary, and can correlate 
the non-neighbouring regions A and B2 directly. Solving 
the problem therefore requires modifying the calculation 
of the topological entropy such that local correlations 
do not cause unsuppressable corrections, whatever their 
range. To do this let us define the topological entropy 7 
as the information that the region A shares with B, but 
not Bi or B2 individually. 
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This definition isolates the correlations that can be 
seen only when all the three regions of the annulus are 
considered. Only correlations that contribute to Ia,b, 
but not to Ia,b ± , Ia,b 2 ° r Ibx,b 2 wn l contribute to 7. 
Since non-topological correlations characterized by two- 
point correlation function will indeed contribute to the 
latter three mutual informations, they cannot affect the 
value of 7, and so cannot cause ambiguity whatever their 
range. Only truly non-local correlations, such as those 
existing along strings in topological ordered states, will 
therefore remain. Note that this definition only requires 
that the three regions together form an annulus, and the 
composite of any two regions does not. As such partition- 
ing the closed loop according to Fig. QJb) is therefore no 
longer a requirement (though it is still valid, and we will 
continue to use in this study), and other tripartions of 
an annulus may be used. 

The only caveat is as follows. Suppose that infinite 
range non-topological correlations are present in a state 
such that the operators Oi are not defined on a single spin 
i, but of spins within its neighbourhood also. Consider 
then a point j that lies on the boundary between the 
regions B\ and B2, such that the support of Oj is split 
between these two regions, and consider also a point i 
that lies fully within the region A. Any two point correla- 
tions between i and j can therefore only be revealed when 
all three regions are considered. Such correlations would 
therefore appear to contribute to Ia,b without contribut- 
ing to Ia.Bi, Ia,b 2 or Ibi,b 2 , and so cause ambiguities 
with the value of 7. However, assuming the Hamiltonian 
causing the correlations is local, the points i and j will 
not be correlated in isolation. There will be other points 
k with which i shares exactly the same information it 
shares with j. In the large annulus limit, at least some of 
these will exist such that that Ok has full support within 
B\ or B-2- The contribution to Ia,b will also be present 
in I a. Bi or Ia,b 2 , respectively, and so be removed from 
the calculation along with all other non-topological con- 
tributions. 

Calculation of the modified topological entropy is not 
trivial, even when only classically correlated states are 
considered. Any calculation must distinguish the infor- 
mation that A shares with both B\ and B2, and that 
which it shares only with the composite region B. In or- 
der to make the distinction, and calculate the modified 
topological entropy, we can use the state of the system, 
p to define an alternative state p' . This is constructed 
such that pABi = p'abx ana ^ Pab 2 — Pab 2 - The correla- 
tions between the regions A and B\ are therefore exactly 
the same for p' as they are for p, as are the correlations 
between the regions A and Bi- We then maximize the 
entropy S(p' AB ) over all states for which these conditions 
hold. Any additional correlations that A has only with 
the composite region B will then not occur in p' . The 
topological entropy is then, 

7 = Ia,b(p) ~ IaAp')- (6) 
This calculation involves both reconstruction of a joint 



state from reduced density operators and a maximiza- 
tion. As such, computation of this quantity may be diffi- 
cult in general. A simple method of calculating 7 in this 
way does exist if the states are classical, i.e. p takes the 
form p = Ea : 6 1 ,6 2 -P( a ' & i> & 2) I a,b 1 ,b 2 ) (a,b 1 ,b 2 , where 
{|a)} is a basis for the states of A, etc. In this case 
the state p' can then be defined using the probabili- 
ties p'(a, 61,62) = p( a )p(bi\a)p(b2\a), the product of the 
marginal distributions for B\ and B2 conditioned on A. 
This maintains the correlations between A and B\, and 
between A and B 2 , but breaks those between A and B 
alone. 

Even without an efficient means to calculate 7 in gen- 
eral, both upper and lower bounds may easily be de- 
termined, applicable in both the quantum and classical 
cases. To calculate these bounds, we consider the quan- 
tity Ia,b — 7, the information A shares with B\ and B2 
individually but not with B alone. This satisfies the 
bounds, 

Ia,b -7 > max^Bi, Ia : b 2 ), 
Ia,b —J< Ia,b x + Ia,b 2 - 

The first of these come from the fact that Ia.b — 7 is at 
least the information A shares with B\ , or that which it 
shares with B2 ■ Taking the maximum of the two hence 
forms a lower bound. The second comes from the fact 
that adding the information A shares with B\ to that 
which it shares with B2 will count all of the information 
in Ia,b — 7, and may even double count some, and so 
forms an upper bound. Using these, and noting that 7 
cannot be less than zero, it follows that, 

7 > min(JA,B - Ia,B! - Ia,b 2 ,0), 

7 < Ia.b - max(/ yliBl , I a ,b 2 )• (7) 

This gives bounds for 7 in terms of the mutual informa- 
tions for bipartitions. Note that these are exactly the 
kinds of quantities used to evaluate Eq. ([3]) and Eq. ([2]) 
above. As such, computation of these bounds is no harder 
than computation of existing quantities. 

To demonstrate that the modified topological en- 
tropy correctly identifies topologically ordered states, the 
bounds may be calculated for the string- net models 0), 
using the relations from [12]. For these, the entropy of a 
region R is given by, 

5b = -jR log D - n R § log Y (8) 

where jn is the number of disconnected boundary curves 
in R and ur is the number of spins along the boundary 
(see [l2j for more details on the specific models consid- 
ered). The two regions A and B2 are not neighbouring, 
and so Jab' = 3 A + 3b 2 and uab 2 = n-A + ns 2 - As such, 
we find that Ia,b 2 — 0; and so the regions are not cor- 
related. The upper and lower bounds of Eq. ([7]) then 
coincide both with each other, and with the definition of 
7 in Eq. ^ . This gives the topological entropy a unique 
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value for these states, equal to the expected result of 
7 = 2 log D. 

For the state p«, I R , R >(pW) = Ir,r>(p (0) ) + 2 log 2. 
As such, since we know that the topological correlations 
result in I A , B (p {0) ) - Ia,bAp (0) ) = 2 log 2, we see that 
the lower bound of Eq.Q vanishes and the upper bound 
gives 2 log 2. For the state p( 2 \ the lower bound vanishes 
and the upper bound diverges with increasing region size. 

The bounds for pW and p^ therefore do not assign 
a unique value to 7 for in either case, nor can they tell 
us whether 7 = or 7 > 0. Nevertheless, the fact that 
they are not equal is a clear signature of long-range non- 
topological correlations. Once their presence has been 
detected, efforts to remove the ambiguity they cause can 
be made. This is in contrast to previous methods, which 
assign a single value and give no clue as to whether this 
stems from topological correlations, or long range non- 
topological correlations. 

IV. CONCLUSIONS 

Here we have shown that existing definitions of the 
topological entropy and topological mutual information 
can give misleading results, since they cannot completely 
distinguish between topological and non-topological cor- 
relations in all cases. As such, a modified definition of 
the topological entropy, 7, has been proposed. This is 
equivalent to the topological entropy of Levin and Wen 
in most cases. However, unlike existing definitions, 7 



is robust against all non-topological correlations, even if 
they effectively have an infinite correlation length. Up- 
per and lower bounds on this modified entropy have been 
determined, which can be computed with the same com- 
plexity required for previous methods. These bounds will 
coincide in most cases and so provide a unique value for 
7. In the the presence of long range non-topological cor- 
relations, the bounds diverge. This serves as a witness to 
the presence of such correlations, as well as limiting the 
ambiguity they can cause. 

It should be noted that problems arising from infinite 
correlation lengths do not occur only for mixed states. 
Pure eigenstates of gapless Hamiltonians will similarly 
cause ambiguities in the calculation of the topological 
entropy. Though a full treatment of gapless states and 
the ambiguities they cause to the topological entropy is 
beyond the scope of this work, it seems likely that the 
modified topological entropy defined here will be useful 
also in the gapless case. 
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